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field dependence of the phonon thermal conductivity is ascribed to the anisotropic 
superconducting gap structure in Sr2Ru04. 
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1. Introduction 

Sr 2 Ru04 first attracted a lot of interest as a superconductor which has the layered 
perovskite structure without copper |T|. Later on it has become clear that despite the 
crystal structural resemblance to the high-T c cuprate superconductor L^-xSr^CuC^, 
the electronic properties of these two systems are quite different J| . The normal state of 
Sr 2 Ru0 4 is well described by Landau Fermi-liquid model ||, and its superconductivity 
is characterized by the spin triplet pairing and spontaneous breaking of the time-reversal 
symmetry § §, g % 

Study of the thermal conductivity can be very useful for probing the quasi-particle 
spectrum in the superconducting state |§. These measurements played an essential role 
in locating the position of the superconducting gap nodes in high-T c cuprates |9|, |1(J 



and in heavy fermion superconductors fLl| , [12|j . Of note, in these systems the thermal 
conductivity is mostly determined by quasi-particle excitations and is electronic in 
origin, as confirmed in the cuprates by measurements of the Righi-Leduc effect (thermal 
Hall effect) ||13 |. Detailed experimental characterization of the electronic thermal 



conductivity stimulated intensive theoretical studies on the subject [[14], [15|, [16], [L7| and 
at present there is a reasonable understanding of the electronic thermal conductivity in 
unconventional superconductors. The thermal conductivity of Sr 2 Ru04 with the heat 



flow along the conducting plane was also shown to be electronic in origin |L8|, [191], and 



it is well described with a model for unconventional superconductors with line-nodes in 



the superconducting gap [19, 20, pT 



It is significant that Sr 2 Ru04 has a quasi-two-dimensional electronic structure 



||, |22j and consequently, all of the electronic transport properties should be strongly 
anisotropic. Therefore it is of interest to extend the measurements to the configuration 
with the heat flow perpendicular to the plane. Here the low electronic conductivity 
provides a rare opportunity to study the •phonon thermal conductivity in unconventional 
superconductors. 

In this article, we report measurements and analysis of the thermal conductivity 
along the least conducting direction (the [001] crystal axis) in the normal, 
superconducting and mixed states of Sr 2 Ru0 4 . Through the analysis of the thermal 
conductivity in magnetic fields and the comparison of the results for the heat flow 
directions parallel and perpendicular to the plane, we separated the phononic and 
electronic contributions and characterized their field and temperature dependences. 

We revealed that the temperature and field dependence of the phonon thermal 
conductivity reflects an increased quasi-particle density in the bulk of unconventional 
superconductors. A sub-linear field dependence of the phonon thermal resistivity is 
ascribed to the presence of nodes in the superconducting gap, which is consistent with 
other relevant experiments; specific heat [23|, NMR [P^] , magnetic penetration depth 



25| and ultrasonic attenuation p6| . 
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2. Experimental methods 

The single crystals of Sr2Ru04 used in this study were grown in an infrared image 
furnace by the floating zone method j27J. In the crystal growth procedure, an excess Ru 
is charged to create the Ru self-flux in order to compensate for Ru metal evaporation 
from the surface of the crystal rod. Therefore, inclusions of Ru metal are frequently 
found in the central part of the crystal boule. Since these inclusions are responsible for 
the formation of the so-called 3K-phase f28|, to obtain the single-phase samples we cut 



them from the surface part of the boule. The absence of the 3K-phase was checked by 
the AC susceptibility measurements. We selected a high quality sample with T c = 1.42 
K, close to the extrapolated impurity - free value T c0 = 1.50 K |29|, [3(J. The sample was 
cut into a rectangular slab elongated in the [001] direction with dimensions 1.5x0.5x0.5 



mm 3 . 



The thermal conductivity was measured by usual one heater - two thermometers 
method with three 2 kQ Ru02 resistance chips ||32j| . The direction of the heat current 



was set along the [001] crystal axis perpendicular to the conducting plane. (Below for 
the thermal conductivity perpendicular to the plane we use the notation k oi, with the 
subscript showing the direction of the heat flow.) The end of the sample was attached 



directly to the cold finger of a miniature vacuum cell |33|. The cell was placed in a 
double-axis sample rotator in a 3 He refrigerator inserted into a 17 Tesla superconducting 
magnet. The rotator enabled precise cell orientation with a relative accuracy of better 
than 0.1°. Simultaneous electrical conductivity measurements were performed using the 
same thermal/electrical contacts. 

3. Results and discussion 

3.1. Temperature dependence of thermal conductivity 

In Fig. D we show the temperature (T) dependence of k oi/T of Sr 2 Ru04 in magnetic 
fields along the [100] direction. k 0Q i/T shows a shoulder on entering the superconducting 
state, contrary to the in-plane thermal conductivity Kioo, of the same high quality (T c 
= 1.44 K) crystals, which shows an immediate decrease below T c |0| . 



To obtain an appropriate understanding of these results, the identification of the 
dominant heat carriers and the scattering mechanisms is important. The thermal 
conductivity k is a sum of contributions of electrons and phonons, k = K e + k, 9 , and in 
the superconducting state both of these quantities change in a complicated way. In the 
normal state, however, throughout the whole temperature range studied we are already 
in the residual resistivity range, therefore, the electronic component n e /T should be 
constant with T. 

The simplest way to estimate K e is to calculate it from the value of the electrical 
resistivity p. The two quantities are related by the Wiedemann-Franz (WF) law, 
px (/t e /T) = L , where L is a Lorenz number of 2.45x 10 _8 Wfi/K 2 . This law was shown 



to be satisfied in Sr 2 Ru04 |19| , contrary to the electron-doped cuprate superconductors 
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[5IJ. To exclude large errors in determination of the sample geometry, we compare the 
measurements of both quantities made on the same electrical/thermal contacts. Using 
the resistance value in the field of 1.5 T along [100] direction, the electronic contribution 
to k oi in the normal state is evaluated as Kq 01 /T = 9.9 mW/K 2 m, which amounts to 
7% of the total Kooi in the normal state at 1.5 K and 32% at 0.3 K. The validity of 
the WF law is based on the hypothesis that the electronic conduction is determined by 
the elastic scattering process. The inelastic processes reduce the thermal current more 
effectively than the charge current, giving the WF ratio smaller than Lq. Thus this 
method of estimation provides the upper bound of the actual electronic contribution. 

Alternative way of estimation of n e uses the assumption that relative changes of 
electrical conductivity p~ l and of n e with fields are the same |34 ]. In this way, using 
the data for magnetothermal resistance and magnetoresistance in the normal state with 
field along [100] direction, we obtain a consistent value Kq QI /T = 9.5 mW/K 2 m, 7% of 
measured k 0Q1 /T at 1.5 K and 30% at 0.3 K. 

As we see above, in the normal state, the thermal conductivity is mainly phononic 
in origin. This is moreover true for the superconducting state, since here, due to the 
dominance of impurity scattering for conduction electrons, the electronic contribution 
n e /T shows notable additional decrease as compared to the normal state fl9|| . Therefore 
we can conclude that the the enhancement of k 0Q i below T c is due to the phonon mean 
free path increase, the cause of which will be discussed below. 

Possible scattering mechanisms in the phonon heat conduction include the 
scattering on the sample boundaries, conduction electrons in the normal state, thermally 
excited quasi-particles and vortex cores in the superconducting state. As can be seen in 
Fig. [I], above T c , k oi/^ changes linearly with temperature (kooi/T oc T), which shows 
that the k 9 is limited mainly by scattering on conduction electrons. Thus the hump 
of Kooi/T observed below T c is interpreted as a result of the condensation of electrons 
into the Cooper pairs which do not contribute to the scattering of phonons. For the 
conventional superconductors, Bardeen, Rickayzen and Tewordt (BRT) calculated the 
increase of k 9 due to the enhancement of the phonon mean free path caused by the 
formation of the Cooper pairs below T c |[35| , assuming an exponential decrease of the 
quasi-particle density. According to their theory, almost 10 times increase of k 9 is 
expected at the temperature of 0.5 T c |]35]| . It can be seen that the actual increase 
of the thermal conductivity is far lower than this expectation, and the large difference 
indicates the existence of additional scatterers such as grain/ sample boundaries or excess 
quasi-particles resulting from the existence of nodes in the gap. 

The role of the boundary scattering can be evaluated in the following way. In 
the boundary scattering regime, the phonon thermal conductivity is described by the 
asymptotic formula k 9 b = ^Cvl oc T 3 , where C=/3T 3 is the phonon specific heat per unit 
volume, v is the mean sound velocity and I is the mean free path of the heat carriers 
limited by the sample dimension. The mean sound velocity along the [001] direction 
of the tetragonal crystal is described by the relation v = Vz ^lt] , where s = vl/vt is 



the ratio of the longitudinal and transverse sound velocities ||36|| . The mean free path 
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of phonons in the boundary scattering regime is Ib = -^=, where d is the mean width 
of the sample cross section |[36| . If we take a phonon specific heat coefficient of (3= 
0.197 mJ/K 4 mole [23], sound velocities vl = 3900 m/s and Vt = 1300 m/s at low 



temperatures |J7| and d = 0.5 mm, we obtain the phonon thermal conductivity in the 
boundary scattering regime as k 9 b = 0.8 T 3 W/Km. This curve is plotted in Fig. |l]with 
a solid line, which lies well above the actual data in the whole temperature range. This 
means that the phonon thermal conductivity is strongly limited by the large amount 
of the thermally excited quasi-particles in the vicinity of nodes. Actually the observed 
height of the hump below T c is 8 times lower than the prediction of the scattering on 
the sample boundaries. 

3.2. Field dependence of thermal conductivity 

Now we discuss the field dependence of the phonon thermal conductivity. For the 
separation of the K e and k 9 , we analyze the different field dependence expected for two 
kinds of heat carriers in comparison with the known field dependence for K e [[TTJ . 

The field dependence of Kqq\/T at 0.3 K in the magnetic field applied along the [001] 
direction is shown in Fig. ^] (a). On field increase from T kqqi/T first remains nearly 
constant, then shows a sharp drop, takes a minimum and starts to increase rapidly 
towards H c2 . Finally Kooi/T flattens above H c2 in the normal state. The plateau 
observed at low fields corresponds to the complete Meissner state with H cl ~ 12 mT at 



0.3 K 38 



Since the effect of magnetic fields on the phonon thermal conductivity comes only 
from the creation of additional scatterers, k 9 decreases monotonically with fields. The 
field dependence of n e of Sr 2 Ru04 is more complicated. At high temperature K e shows 
initial decrease with fields, and then increases towards H c2 [0. However, at low 



temperatures n e increases monotonically with fields |T{|. Thus the rapid increase of 



Koqi/T near H c2 comes from the increase of n e on the slowly decreasing background 
of k? . This allows us to decompose the field dependence into electronic and phononic 
components. We assume that the electronic part in k 001 has the same field dependence, 
except for difference in amplitude, with as /tioo which is purely electronic in origin. Then 
the field dependence of 

^001 Cclll be determined by subtracting the scaled /tioo (i.e. Kqoi) 
from the total Kooi as shown in Fig. RL Here the scaling factor is used as adjustable 
parameter, chosen so as to give coincidence of the scaled Kioo and Kooi curves near 
H c2 . As shown in the figure, in the superconducting state, the phonon contribution is 
strongly dominant, accounting for 88% of Kqqx/T at T and 0.3 K. Although being 
efficient, this procedure is, strictly speaking, not correct because we neglected the 
difference in magnetoresistance and hence magnetothermal resistance between k\ qq and 
Kooi- However, since the magnetic field necessary to suppress the superconducting state 
in H || [001] configuration is small, this correction is not of large importance here. 

On the contrary, in the magnetic fields parallel to the plane, the magnetoresistance 



becomes large for the current perpendicular to the plane even just above H c2 40 . In our 
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case the resistance at H c2 in perpendicular field is approximately 1.7 times smaller than 
in parallel field. This magnetoresistance accounts for the smaller normal state electronic 
thermal conductivity in parallel fields (30-32 % of the total conductivity at 0.3 K, see 
Section |3.1| ) as compared to perpendicular field case, in which case n e constitutes about 
a half of the total normal state conductivity (Fig. |3|). 

The dependence of k 001 /T on magnetic fields H || [100] (along the plane) at 0.3 K 
is shown in Fig. || (b). Contrary to the perpendicular field case, k 001 /T curve does not 
have the minimum, resulting from increase of the electronic conductivity towards H c2 . 
Instead, the conductivity decreases gradually and then shows a rapid decrease towards 
the normal state value in a narrow field range in the vicinity of H c2 . We point out here, 
that the absence of the upturn at H c2 in parallel fields at low temperatures, despite 
a notable increase of n e [fL9|l , is completely unusual. This fact implies that (i) phonon 
scattering is strongly suppressed in the superconducting state at low fields, (ii) some 
additional mechanism for phonon scattering appears close to H c2 , which overcomes a 
rapid increase of K e . The latter fact cannot find any explanation in the standard vortex 
scattering scenario. As is clear from the comparison with Fig. |3], the absence of the 
upturn rules out the application of the procedure we used in the perpendicular field for 
extracting n 9 (H) in the superconducting state. 

As we already mentioned above, the field dependence of the phonon thermal 
conductivity is totally determined by the variation of the scattering rate. The scattering 
in the superconductor is governed by the vortex cores and delocalized quasi-particles 
in the bulk, the latter acting similar to conduction electrons in the normal state. For 
a superconductor with an anisotropic gap with line nodes the quasi-particles in the 
bulk are generated by a magnetic field, as was first pointed out by Volovik |^TJ. This 
effect is caused by a supercurrent flow around a vortex and associated local Doppler 
shift of the quasi-particle spectrum. The shift causes influx of the quasi-particles from 
the core into the bulk along the directions of the nodes in the superconducting gap. 
This increase of the quasi-particle density in fields is indeed observed in specific heat 
measurements of Sr 2 Ru04 p3| , giving a strong evidence for the nodal structure of the 
superconducting gap. At low temperatures the density of delocalized QP is negligible 
in conventional type-II superconductors with the isotropic gap, therefore the phonon 
scattering in magnetic fields between H c \ and H c2 proceeds predominantly on the vortex 
cores. In this case the phonon thermal resistivity W(H) increases linearly with H for 
H > H c i, being proportional to the number of vortices. It can be expressed |fl2" 



as 



W(H) = — = W + W n H/H c2 , (1) 

where Wo and Wq + W n represent the phonon thermal resistivity in the Meissner and 
normal states, respectively. 

Fig. f| (a) shows the field dependence of the phonon thermal resistivity denoted as 
W 9 (H) in the magnetic field H || [001]. The phonon part was extracted as shown in 
Fig. |3|. W 9 (H) increases linearly with field and shows the dominance of the vortex core 
scattering. In this field configuration the coherence length is large (£ a & ~ 660 A)p9[], 
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giving a large cross-section of the vortex core for the phonon scattering even at low 
fields. Thus the delocalized quasi-particles, despite their density is large as seen in the 
specific heat measurement [23], do not effectively influence the phonon transport. 

In Fig. [| (b) we plot the thermal resistivity W(H) at 0.3 K as a function of magnetic 
fields parallel to the plane (H || [100]). Here W(H) is obtained directly from l//cooi and 
consequently contains a small electronic contribution. However, for fields H < H c2 /2, 
the electronic contribution is slightly varying with fields ||19| , and should give small and 
smooth offset of zero line, hence the shape of W(H) can well represent that of W 9 (H). 
As it can be seen from the figure, the W(H) shows a notable sublinearity and at low 
fields exceeds the if-linear increase expected for scattering on the vortex lattice. This 
shows that the vortex scattering is not dominant in this configuration, which can be 
understood if we recall that the vortex core for parallel field is approximately 20 times 
smaller (£ c ~ 33 A) [39 than in the perpendicular field. In the situation when the vortex 
scattering is weak, the delocalized quasi-particles induced by magnetic fields can become 
the main scatterers of phonons. 

A quantitative description of the role of the delocalized quasi-particles on the 
phonon transport is very difficult at present since it is almost unexplored theoretically. 
As a naive discussion, however, we can assume that the thermal resistivity is proportional 
to the number of scatterers, i.e. the delocalized quasi-particles. This density was 
calculated by Volovik for isolated single vortex core in the <i-wave superconductor [ffl}| 
and it is proportional to \/H at very low temperatures. Thus we may consider that 
the phonon resistivity should be proportional to \/~H . The curve in Fig. [4] (b) is a 
most close fit of the data with a v^ff-function. While this fit reproduces our result well, 
the above discussion is oversimplified and a proper theory for the phonon transport 
in the unconventional superconductor is strongly desired. Very recently, Won and 
Maki ||3| calculated the field dependence of the phonon thermal resistivity in the nodal 
superconductor, assuming dominance of the longitudinal phonons in the heat transport 
at low temperatures and the validity of the Matthiesen rule for rates of scattering on 
sample boundaries and on quasi-particles. They derived a sub-linear field dependence 
of the form HlnH due to the delocalized quasi-particles |[43|| , which is not very distinct 
from sqrt(H) dependence at low fields and seems consistent with our results within the 
experimental error. 

To check the relation between the phonon resistivity and the quasi-particle density, 
we compare in Fig. [5] the field dependence of the thermal resistivity and of the specific 
heat by Nishizaki et al. representing the quasi-particle density, both taken at 0.3 K 
in the same experimental configuration with in-plane fields. Note that the horizontal axis 
is scaled as y/H. It can be seen that at low fields both quantities show similar sublinear 
dependence and match well, supporting our assumption on the phonon scattering on 
the delocalized quasi-particles. Despite this similarity, the two curves can never coincide 
completely, since corrections due to the electronic contribution and additional scattering 
on the vortex lattice are essential for the thermal conductivity at high fields. When 
scaled to overlap in the low field region, the thermal resistivity appears to be higher than 
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the specific heat curve. The electronic contribution n e , which can become important 
near H c2 , should cause the decrease of the resistivity and thus can not be the reason 
for this deviation. Therefore we conclude that additional scattering on vortex lattice is 
important as well. 

An additional support for the role of scattering on field-induced quasi-particles in 
Sr 2 Ru0 4 in parallel field comes from an unusual shape of the k(H) curve near H c2 at low 
temperatures. As can be deduced from Eq. (1), in fields close to H c2 the decrease of the 
phonon contribution due to scattering on the vortex lattice saturates to a value equal 
to that in the normal state. Therefore, as in the perpendicular field case, the increase 
of the K e towards H c2 over this flat background gives an increase of total k. However, 
the situation in Sr 2 Ru04 in magnetic fields parallel to the plane is notably different. 
Due to a multi-component order parameter, the superconducting state of Sr 2 Ru04 
is characterized by existence of at least two superconducting phases [f44| . The main 
superconducting state occupies almost all H — T domain of the superconductivity, while 
the second superconducting state exists at low temperatures (below 0.8 K) at precisely 
aligned (within ±3° close to the plane) fields in close proximity to H c2 |Tj| [23[ [45|j . 

A phase transition from low field superconducting phase to high-field one is 
accompanied by a release of large density of quasi-particles. This release of QP gives a 
rapid increase of electronic thermal conductivity, which in the normal situation should 
give an upturn in k,(H). However, since the release of the quasi-particles is similarly 
important for the phonon scattering, and the phonon contribution k 9 is much larger 
than K e , the decrease of k 9 is dominating over increase of n e and gives downturn in 
Koai(H). 

To see whether the unusual curvature is indeed related to the formation of the 
second phase we show in the inset of Fig. ^ the k(H) dependence at 0.82 K, 
above the domain of the second superconducting phase. As can be seen, the curve 
shows usual shape with a small upturn in k(H) near H c2 , despite notably smaller 
electronic contribution to the total thermal conductivity at this elevated temperature. 
Simultaneously, as shown in the main panel of Fig. ||, the downward curvature appears 
in the same range of field inclinations (within ±3° to the plane) where the second 
superconducting state exists and where rapid release of quasi-particles is observed 
in specific heat and electronic thermal conductivity measurements [[l^, These 
observations clearly show that the downturn is determined by formation of a second 
superconducting phase, and is determined by a rapid increase of scattering on quasi- 
particles near H c2 . 



4. Conclusion 



The thermal conductivity perpendicular to the conducting plane of Sr 2 Ru0 4 down to the 
lowest temperatures is dominated by phonon contribution, with scattering by conduction 
electrons or quasi-particles as the dominant scattering mechanism. The shoulder in the 
temperature dependence of the thermal conductivity in the superconducting state is 
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related to the enhancement of the phonon mean free path due to the condensation of 
normal electrons below the superconducting gap, and its height is substantially lower 
than the expectation for a superconductor with a nodeless gap. In the mixed state 
the phonon heat transport is determined by the scattering on the vortex cores and on 
the bulk quasi-particles outside the core. The latter contribution is important in the 
configuration when the magnetic field is parallel to the plane, where the core radius is 
small. These results imply that the phonon thermal conductivity is a useful probe for 
the characterization of the unconventional superconducting state, complementary to the 
electronic thermal conductivity or specific heat E6| . 
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Figure 1. Temperature dependence of the thermal conductivity perpendicular to 
the highly conducting plane, kqoi, of Sr2Ru04 in magnetic fields of T, 0.3 T and 1.5 
T (normal state) along [100] direction. A solid curve is the calculated phonon thermal 
conductivity in the boundary scattering mode. 
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Figure 2. The magnetic field dependence of k 00 i/T at 0.30 K; (a) in magnetic fields 
along the [001] direction, (b) in magnetic fields along the [100] direction. 
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Figure 3. The magnetic field dependence of kqqi/T in magnetic fields along the [001] 
direction at 0.30 K (closed symbols). Open symbols represent the field dependence of 
the thermal conductivity with the heat flow along the conducting plane Kioo, measured 
in the same field orientation and temperature, scaled to fit the increase of Kooi near 
H C 2 . The solid curve shows the difference, representing the field dependence of the 
phonon conductivity Kq 01 . 
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Figure 4. (a) The dependence of the phonon thermal resistivity on magnetic fields 
perpendicular to the plane (H \\ [001]), extracted in the way shown in Fig. ||. (b) The 
low field part of the field dependence of the thermal resistivity in magnetic fields along 
[100] direction within the plane (symbols). The line represents a least square fit of the 
data with a \/i7-function. 
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Figure 5. The field dependence of the thermal resistivity W(H) (closed symbols) and 
specific heat (open symbols) in Ref. [23] in magnetic fields along [110] direction within 
the plane at 0.30 K. The horizontal axis is scaled as \Jli. 




Figure 6. Transformation of the field dependence of kqqi/T at base temperature 
(0.30 K) with direction of the field. The numbers show an angle of field inclination to 
the conducting plane, being 90° for parallel alignment; the [110] direction. Inset; the 
magnetic field dependence of kqqi/T at several temperatures with H applied parallel 
to the conducting plane ([110]). Note that some offset is subtracted by each curve. 
Arrows show H c2 for each curve. 



